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Abstract 

Flow forming involves complicated tooling/workpiece interactions. Purely 
analytical models of the tool contact area are difficult to formulate, resulting 
in numerical approaches that are case-specific. Provided are the details of an 
analytical model that describes the steady-state tooling/workpiece contact 
area allowing for easy modification of the dominant geometric variables. The 
assumptions made in formulating this analytical model are validated with ex- 
perimental results attained from physical modelling. The analysis procedure 
can be extended to other rotary forming operations such as metal spinning, 
shear forming, thread rolling and crankshaft fillet rolling. 

Keywords: flow forming, metal forming, physical modelling, contact 
interface, analytical model 



Nomenclature 

a - the entry angle from workpiece to the roller, also known as the attack 
angle [°] (Fig. 4) 



A xyz - overall contact area [mm 2 ] (Section H]) 
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Axy - %y planar projection of the contact area [mm 2 ] (Section H]) 

A xz - xz planar projection of the contact area [mm 2 ] (Section @| 

A yz - yz planar projection of the contact area [mm 2 ] (Section H]) 

(3 - the trailing angle from the workpiece to the roller, also known as exit 
angle or planishing angle [°] (Fig. 4) 

5d - distance between analytical surface and experimental surface nearest 
neighbour points (Section [3]) 

Sr>i - distance between analytical surface and experimental surface inter- 
polant (Section [3]) 

d - the distance between the center of the mandrel and the center of the 
roller [mm] (Fig. 4, d = R m + t f + R + R r ) 

f z - axial feed rate of the roller down the face of the cylinder, along the z 
direction [mm/min] 

MSE - Mean Square Error (Eq. El]) 

n - the mandrel rate of rotation [revolutions/min] 

P - the roller path pitch or distance traveled axially by the roller in one 
revolution [mm] (P = f z /n) 

R - roller nose radius [mm] (Fig. 4) 

R* - numeric resolution of the solution (Section I2.2.2P 

Ri - initial workpiece radius (R4 = R m + t ) 

R m - the mandrel radius [mm] (Fig. 4) 

R r - the roller radius excluding the radius of the nose [mm] (Fig. 4) 
S - Intermediate set of radial quantities used to find X p and Y p (Eq. 126]) 
6f - the angle of contact between the roller and workpiece [rad](Fig. 2, Eq. 



2 



6i - intermediate value of 6 used for the iterative solution of the contact area 
[rad] (Eq. [33]) 

6 max - maximum angular limit of the solution (Eq. [T4"]) 

@ixR* - angular coordinates used to define boundary surfaces (Eq. [T9]) 

to - starting material thickness [mm] (Fig. 4) 

tf - the final material thickness [mm] (Fig. 4) 

Xi - x coordinates lying on the instantaneous roller position (Eq. [201) 

Xi - x coordinate used for intersection conditioning (Eq. [5]) 

x\ - x coordinate used for intersection conditioning (Eq. [3J) 

X m - x coordinate lying on the cylinder defined by Ri (Eq. 123]) 

i max - maximum limit in the x direction of the solution (Eq. [T5]) 

X p - x coordinates lying on the previous roller path (Eq. K7j) 

X s - x coordinates within the roller /workpiece contact area (Eq. [30]) 

x u - x coordinate used for intersection conditioning (Eq. [lj 

Yi - y coordinates lying on the instantaneous roller position (Eq. |2T]) 

Y m - y coordinates lying on the cylinder defined by Ri (Eq. [24]) 

?/max - maximum limit in the y direction of the solution (Eq. [T6]) 

Y p - y coordinates lying on the previous roller path (Eq. [28]) 

Y\xr* - y coordinates used to define boundary surfaces (Eq. ITS]) 

Y s - y coordinates within the roller /workpiece contact area (Eq. [31]) 

Zi_ 2 - axial limits of the workpiece/roller contact area, z coordinate of the 
endpoint of contour 1 and starting point of contour 2 (Eq. [T2land[T3]) 

Zi_3 - axial limits of the workpiece/roller contact area, z coordinate of the 
endpoint of contour 1 and starting point of contour 3 (Eq. [71 to [TO]) 
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Zi - z coordinates lying on the instantaneous roller position (Eq. 
Zi - z coordinate used for intersection conditioning (Eq. [6]) 
Zi - z coordinate used for intersection conditioning (Eq. H]) 
Z m - z coordinates lying on the cylinder defined by Ri (Eq. 1251) 
Z p - z coordinates lying on the previous roller path (Eq. 129]) 
ZixR* - z coordinates used to define boundary surfaces (Eq. ITTj) 
Z s - z coordinates within the roller /workpiece contact area (Eq. I3"2"j) 
z u - z coordinate used for intersection conditioning (Eq. [2]) 

1. Introduction 

To determine the energy required to form a component, the size and 
orientation of the tooling interface on the workpiece is necessary. While 
purely analytical models describing this contact are preferable, they are usu- 
ally difficult to attain for complex metal forming processes. In this study, 
an analytical approach is presented to model the tooling/workpiece contact 
area in an application of rotary forming. While the present work focuses on 
an implementation for flow forming, the applied technique can be applied 
to other variants of rotary forming operations such as metal spinning, shear 
forming, thread rolling and crankshaft fillet rolling. 

Flow forming, a variant of metal spinning, is a process us ed to fabri 



cate r otationally symmetrical parts from ductile materials, after IWong et al 



(120031 ). During flow forming, the workpiece is clamped to a rotating mandrel 
and pressed into contact with the mandrel by rollers. The rollers induce high 
levels of plasticity in the workpiece causing it to undergo both reduction in 
thickness and axial lengthening. Since the rollers press on only a very small 
area of the overall workpiece at any given time, the deformation is highly 
localized between the roller and workpiece. To properly understand the dis- 
tribution of this intense local plastic deformation it is essential to be able to 
calculate the roller /workpiece contact area from the geometric parameters 
that govern the flow forming process. In addition, the roller /workpiece con- 
tact area is critical to coupling other experimental findings, such as power 
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consumption, frictional effects, force, stress and strain distributions through 
the workpiece back to geometric process parameters. 

In flow forming, the combined mandrel rotation and linear movement 
of the rollers induce contact on the workpiece along a helical path. This 
helical tool path, coupled with the curved profile of the rollers leads to a 
very complicated roller /workpiece contact area. 

In terms of related tool contact studies, an important analy tical deriva 



tionof the workpiece contact in shear spinning was completed by lChen et al. 



(200 51). However, in a comprehensive review of metal spinning processes, 



Music et al.l ( 120101 ) highlighted that the mechanics of flow forming are quite 
different than shear spinning. This is also true for the contact area formula- 
tion as there is little roller penetration into the workpiece and deformation 
proceeds according to the sine rule. 

In terms of flow forming specific research, investigations made by Gur and 
Tirosh (1982), Singhal et al. (1995), Ma (1993) and Jahazi and Ebrahimi 

(2000) have proposed analytical models of this contact. Gur and Tirosh 
(1982) developed the formulation of a planar contact area in each of the pri- 
mary rolling and extrusion deformation directions in backwards flow forming. 
Singhal et al. (1995) derived the contact area imposed by tooling in the flow 
forming of small diameter tubes where the assumption made was that mate- 
rial is assumed to be perfectly plastic, and the tools were assumed rigid. Ma 
(1993) extended the work of Gur and Tirosh (1982) to derive a critical angle 
of attack and Jahazi and Ebrahimi (2000) extended the contact formulation 
made by Gur and Tirosh (1982) to investigate the mechanics in a specific 
application of flow forming. More recently, Kemin et al. (1996), Xu et al. 

(2001) and Hua et al. (2005) have developed Finite Element (FE) models 
of single roller flow forming. In each of these studies, contact was modeled 
explicitly within each respective FE m odel. Furthermor e, with the exception 



of the work by IXu et al.l fcOQlh and iHua et all (l2005h . all previous works 



have made assumptions concerning the roller /workpiece contact geometry 
that do not necessarily reflect the actual contact during flow forming. These 
assumptions include: 

1. Idealized roller geometry (i.e. no blending radii) (Fig. 1); 

2. The use of two-dimensional treatments that do not account for the 
three-dimensional aspects of the workpiece contact; 

3. Not considering the influence of prior forming steps (i.e. roller path 
overlap) on the instantaneous roller /workpiece contact area. 
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The most successful technique for modelling the roller /workpiece contact 
area, and other fa c ets o f the flow forming process, has been through FE 
analyses. Ku et al.l (j200lf ) addressed items 1 and 2 listed above in their work 
to numer ically calculate t he roller /workpiece of a single roller flow forming. 
However, IXu et al.l (j200l[ ) did not give the details of t heir calculation o f the 
contact area, nor did they specifically address item 3. iHua et al.l ( 120051 ) has 
developed a thorough 3-D FE model that addresses all three items, but an FE 
approach is still limited to case-by-case application involving extensive pre- 
processing and explicit geometric modelling. An analytical solution provides 
a solution with significantly lower effort. In the present work, a generalized 
solution is developed for the roller/workpiece contact area during a single 
roller flow forming operation that accommodates items 1 to 3. To accomplish 
this, the following assumptions are made: 



1. 



2. 



The single roller flow forming process proceeds under steady state con- 
ditions. The final and starting thickness, mandrel rotation and feed 
rate are constant. 

The deformation response of the workpiece is perfectly plastic. Elastic 
effects are not considered. 

Volume of the flow formed workpiece is conserved outside the tool in- 
terface. 

No material build-up occurs in front of the roller as the workpiece 
conforms completely with the rigid roller. 



2. Contact Solution 

During flow forming, the roller contacts the workpiece along a path having 
a constant pitch (Fig. 1). The profile of the roller can be divided into three 
regions; the entry region, the nose region and the exit region. These regions 
dictate the size and shape of the roller/workpiece contact area. The contact 
area is bounded by three contours: the tangential exit contour, the axial 
entry contour and the axial exit contour, labeled 1-3 respectively in Fig. 
2(a). The contact area extends angularly from the tangential exit contour 
(6 = 0) through to 9 = Of (Fig. 2(b)). If the roller has an archetypal 
flow forming profile similar to that shown in Fig. 1 with distinct flat entry 
and exit regions and a blending radius between the two that creates a nosed 
roller, the final contact area is dependent on six surfaces (Fig. 3). Contour 
1, and the starting points of contours 2 and 3 (Fig. 2) can be calculated 
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Figure 1: Single roller contact in flow forming showing the mandrel and key roller profiles. 
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Figure 3: Roller profiles deciding the instantaneous contact area during flow forming. 
There are six in total: the nosed region of the roller from the previous workpiece rotation, 
the entry region of the roller from the previous workpiece rotation, the instantaneous roller 
exit region, the instantaneous roller nosed region, the instantaneous roller entry region and 
the outer surface of the unformed workpiece. 

directly as they lie exclusively on the xz-plane. The a priori z-coordinates 
of the extents of contour 1 define the axial limits of the roller /workpiece 
contact area. Once the a solution has been found for the starting and ending 
points of contour 1 (by definition the starting points of contours 2 and 3), 
the common end point of contours 2 and 3 is then solved using an iterative 
technique. 

2.1. Axial Limits 

It is first necessary to calculate the axial limits of contact by determining 
the endpoints of contour 1. Contour 2 is a function of the instantaneous roller 
contact with the workpiece at pitch P = 0. Contour 3 is a function of the 
instantaneous roller contact on the material as well as the tool contact on the 
workpiece one revolution of the mandrel beforehand, at P = f z /n. Contour 
1 exists solely on the xz plane and is bound by the points of intersection 
with contours 2 and 3. Contour 1 is both dependent on roller geometry 
and the roller path pitch, P. There are four possible conditions describing 
the intersection of the current roller position with that of its position on 
the previous mandrel revolution (Fig. 4). Calculation of the location of 
the upper end point of contour 1 for the four conditions shown in Fig. 4 is 
accomplished through comparison of the endpoints of the roller nose profile 
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Figure 4: Upper limit conditions of the contact area. Position P = corresponds to 
the instantaneous roller position and position P = f z /n corresponds to the roller at one 
mandrel revolution beforehand. The upper endpoint of contour 1 can occur within the 
nosed region of the roller on both the instantaneous position and the position on the 
previous mandrel revolution (condition A). It can also occur at the intersection of the 
exit/entry profiles (condition B), the nosed/entry profiles (condition C) or the exit/nosed 
profiles (condition D) of the instantaneous and the previous roller positions. 

on the xz plane. For comparison purposes, the local coordinate system is 
moved on the x axis from the global origin by d — R r — R (Fig. 5). The x 
and z coordinates of the upper end point of the nosed region of contour 1, 
x u and z u (Fig. 5): 

x u = R(l - cos/3) (1) 

z u = Rsin/3 (2) 

For the lower x and z coordinates of the end point of the nose region of 
contour 1, xi and Zi (Fig. 5): 

xi — R (1 — cos a) (3) 

z\ = — .Rsina (4) 

The entry profile of the previous roller path and the instantaneous exit profile 
of the roller will occur at Xi and Z{. These are expressed as: 
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Figure 5: Two dimension cutaway of a single flow forming operation showing critical 
geometric variables and forming zones. 



-Rfsin a cos (3 + sin (3 cos a — sin (3 — sin a) + P sin a sin (3 ,_. 

Xi = T~- i • a 

cos (3 sin a + sin (3 cos a 



R(cosa — cos/3) + P cos /3 sin a 

2j = (6) 

cos p sin a + sin p cos a 

The values of x u , xi, Xi, z u , z\ and Z{ can be compared to identify which con- 
tact condition shown in Fig. 4 applies. The conditions and the relationships 
that must be simultaneously satisfied are shown in Table [TJ 
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Tabic 1: Upper axial limits of contact 



Condition 


Relationships 


Intersection 

(P = 0/P = f z /n) 




Zi + P<Zi 




A 


%u %i 
Xi > Xi 
X u ^> Xi 


Nosed region/Nosed region 




Z t + P> Zi 




B 


Xi < Xi 
X u < Xi 


Exit region/Entry region 




Zl + P>Z t 




C 


%u ^* %i 
X\ < Xi 
X u ^> Xi 


Nosed region/Entry region 




Z t + P<Zi 




D 


%u %i 
X\ > Xi 
X u ^ Xi 


Exit region/Nosed region 



Once the proper contact condition is determined, a solution for the z 
coordinate for the upper end point of contour 1, as well as the upper axial 
limit of the solution space, -21-3, is possible. Solving ^i_ 3 for the appropriate 
condition A through D: 



*i-3 {A} = \ (7) 
Zl - 3 {B} = Zi (8) 



_ r z = v^ 2 -(x-r) 

^ tan a V tan a 



2 

Zl + P) 



(9) 



-3{D} = {^i^i^} (10) 

I z tan/3 ^ u tan/3 J 

In conditions C and D, Z1-3 is expressed as the solution that satisfies the two 
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equations for z. 

The lower end point of contour 1, Zi_2, occurs at the intersection of the 
profile of the instantaneous roller position and the cylinder with radius Ri. 
This intersection depends on roller geometry and the depth that it penetrates 
into the workpiece. Either the roller intersects the workpiece at the flat entry 
region (condition I) or the nosed region (condition II). If the roller intersects 
at the flat entry region: 



where xi is given by Eq. [3j Otherwise, condition II prevails and the roller 
intersects within the nosed region. The solution for this intersection point 
yields the following expressions for the lower end point of contour 1, and the 
lower axial limit of the solution space. This value, Zi_2, for each condition is 
given as: 



2.2. Solution Boundaries 

Thus far, the region where the tool contact resides has been explicitly 
bound in the axial direction between zi_ 3 and £1-2 ■ The following describes 
how the components needed for a computation of the full three dimensional 
contact area are developed. These components are the maximum angular 
limit that the solution space can be defined by and the surfaces that bind 
the solution space of the instantaneous contact. 

2.2.1. Maximum Angular Limit 

Contours corresponding to the ones described in Section [2] that pass 
through the axial limits, ^i_3 and zi_ 2l are formulated to extend angularly 
from 9 = to 9 = 9 max . This value is the absolute maximum value that 9f 
can be, corresponding to P ~ 0. The extremal point at angle (9 max lies on 
the xy plane at z = 0; its coordinates are obtained using the same derivation 
as the general solution for the contact of two circles using the global datum. 
The first circle is one centered at x, y = with a radius of Ri and the other 
is at a distance x = d and y — with a radius of R r + R. 



Xi < t f 



(11) 




(12) 




(13) 



12 



#ma X = arctan ( ) (14) 



with x max and y mSLX as follows: 



_ d^-{R T + Rf + R} 

■^max — 2^ \ ) 



A<PR*-{<P-{R r + RY + Riy 

2/max = ^ (16) 

2.2.2. Surface Definitions 

As described in Section |2j there are six different surfaces that define the 
area. Due to the system complexity, a numerical technique is employed to 
generate these boundary surfaces in three dimensions. In this technique, the 
overall solution space is represented by a finite number of points or nodes. 
For example, if the solution space was broken up into 20 x 20 x 20 uniformly 
spaced points, then the resolution, R*, of the space would be 20. Figure 6 
shows the effect on the solution by increasing or decreasing R* . The following 
are the definitions of the arrays of coordinates used to define the boundary 
surfaces as functions of R*, where initially n = R* and decays for every term 
included in the array until n — 1. 



Z lxR * (n) = z u - (R* - n) [~^—{ J for n = R\ R* - 1, . . . 1 (17) 

>W (n) = (R* - n) for n = R\ R* - 1, . . . 1 (18) 

6i xR * (n) = (R* - n) — ^— for n = R\ R* - 1, . . . 1 (19) 
K — 1 

Now that the arrays of points are formulated, the boundary surfaces can 
be formed as m = R* by n = R* square matrices for each direction through 
space to form Z , Y and 9 . The m direction of these matrices is a 

i m,n ' m,n m,n 

solution for z at a given value in the n direction of y or 9, corresponding 
to the coordinate arrays given in Equations [T71 [18] and [191 The following 
are functions of discrete entries in the matrices that define the boundary 
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(a) 



(b) 




y 



(c) 

Figure 6: Graphical progression of the iterative algorithm used to solve the contact area. 
Starting with a coarse result shown in (a) with R* = 10, with and resolution increased in 
(b) to R* = 40, and finally with a high resolution answer in (c) with R* = 160. 
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surfaces. 

Starting with the instantaneous roller position: 



X, 



^ , i?(cos 0— cos a) / 
cos a cos P v 

d-{<i>-Z m ^(i)J\- *% B ™ 



Z tan ft) a/ 1 - ^ si "7 m ' n Z 

m,n " / V 0^ m,n u 



R 2 -Z 2 



z-,, > Z > z/ 

U — m,n — t 

(20) 



where = R r H — — . Zj and K remain the matrices of dimension R* 

' COS CK L m,n <>m,n 



by iT 



Y —Y 

1m, n m 

Z; _ = Z_ 



(21) 
(22) 



The cylindrical surface defined by Ri that describes the outer surface of the 
workpiece formed around the mandrel: 



X m = Rj cos 9 


(23) 


Y m = Rj sin 9 

i"'m,n t- m,n 


(24) 




(25) 



For the previous roller path, the surfaces are most easily defined in radial 
expressions, and translated to cartesian coordinates. This radial quantity, S 
is defined as: 



' ^ _|_ R(cos f}— cos a) 

cos a cos P 



s= < 



• • • (^tan/3 - Z_ + P (l - %f))) 0- 1 



d + tan a 



tan p 



+ P[1 



2?r 



k d-Rr-y/^-(z m . n -P(l- 9 -^))' 



Z m , n >z u + P 
Z m , n < Zl + P 

Z U + P> ^ >Z t + P 



(26) 

where = R r + Converting these radial quantities into cartesian coor- 
dinates results in: 

(27) 
(28) 
(29) 



X P m, n =Scos9 n 



Pm,n 



z = z = z 

Pm,n Ttlm,n &m,l 
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2.2.3. Contact Surface Solution 

Once the boundary surfaces have been defined, the contact surface can 
be solved. The contact surface points reside within the instantaneous roller 
position definition, bound by the intersections with the workpiece surface 
and the previous roller path surface. In order to determine which points lie 
within these boundaries, the matrices containing the surfaces are parsed with 
logical operators and placed in a new set of matrices corresponding to the x, 
y and z points lying on the contact surface. 

For the x-coordinates of the contact surface, X s belongs to the set of 
coordinates corresponding to: 

X s = XAXi <X m AXi <X V } (30) 

Similarly, for the y-coordinates, Y s : 

Y s = YAYi <Y m WYi <Y V } (31) 

where 'A' is the conjunction (and) operator, and 'V' is the disjunction (or) 
operator. Finally, for the z coordinates, Z s : 

Z, = Z i {Z imin <Z Pm J (32) 

Each of these conditions must be simultaneously satisfied for X s , Y s and 
Z s . While all x-coordinates will lie between the cylinder defined by Ri and 
the previous tool path (leading to the 'A' operator), due to the curvature 
of the instantaneous roller position versus the counterclockwise curvature 
of the workpiece, the ^-coordinates require the 'V' operator for the same 
comparison. 

2.3. Iterative Calculation of the Contact Length 

The contact area extends from 9 = on the xz plane to some final radial 
distance at 6 = 9f. Depending on the geometry involved, the resolution 
selected and the initial value for 9 selected, it is possible that only half of the 
coordinates in X it Y^ and Z^ correspond to X s , Y s and Z s . To overcome this, 
it is necessary to iterate on the initial value of 9 defining the solution space via 
the following method. The initial solution without any iteration (i = 0) has 
all of the surfaces defined extending axially through the endpoints £i_3 and 
Z1-2 of contour 1 (Fig. 2) and angularly from 9 = to 9 i=0 , corresponding 
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according to #j = o = # max - Once the first values for X s , Y s and Z s are available, 
a preliminary calculation for Or is possible, residing within Y s as per Eq. 



9ftt9 i = s > m\ — — — I (33) 

After priming, 6^=0 is probably close to the actual value of Of. However, 
the solution space that defines 9 i= i is spread uniformly across the range of 
< 9 < 9 max and therefore the solution should be retried at a smaller value 
than 9 = 9 max . This can be continued n times until 9i =n is a large percentage 
(~ 95%) of 9 i=n _i to ensure the most accurate result at a given value of R*. 

Once the coordinates for the solution surface have been solved, the total 
surface area can be readily calculated. This can be accomplished through 
tessellating or meshing the coordinates within X s , Y s and Z s , then using 
Gaussian quadrature or a brute force method to find the area of each element 
and summing. 

3. Experimental Validation 

Flow forming is typically applied to forming metallic components involv- 
ing high speeds, feed rates and forming forces. Due to the nature of the 
process, it is difficult to completely stop the process at a particular point in 
time in a safe manner such that the instantaneous forming zone is preserved. 
Therefore, the forming speeds, forces and feed rates must be significantly 
reduced in order to study the tooling/workpiece interaction in flow forming. 

In order to validate the analytical model above, a physical model of flow 
forming was developed using Plasticine conforming to ASTM D-4236. The 
suitability of using Plasticine and simil ar compounds to mo d el me tal form- 



in g processes has been estab lished by ISofuoglu and Rastyl (120001 ) as well 



as iPertence and CetlinI (1997 ). All material preparation steps detailed by 
Sofuoriu and Rastvl (I200OT were followed. The contact condition used for val- 



idation was designed to provide geometry that other modelling efforts have 
failed to address. Specifically, the contact of a nosed roller with a nose/nose 
entry condition of type 'All' (Section I2.ip . Unlike the other contact con- 
ditions, this type is the most complex in terms of curvature as there is no 
straight/linear section appearing anywhere on the surface. 

The tooling used was a smooth mandrel with R m = 69.33 mm and a 
roller with R r = 56.23 mm, R = 5.00 mm, a = 45° and = 60°. These 
were installed on a lathe with a thread-cutting feed set such that P = 2.54 
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mm/rev. The mandrel was dusted with talcum powder and Plasticine of 
a uniform thickness t a = 8.05 mm was set on it. The outer surface of the 
Plasticine was also dusted with talcum powder. The roller was brought into 
contact with the Plasticine such that tf = 6.17 with n = 5 rev/min. The steel 
and Plasticine process components were all measured with standard contact 
measurement apparatus with an accuracy of ±0.01 mm. This reduction level 
was selected such that it maximized the size of the contact patch with the 
given tooling and minimized bulging of material ahead of the roller. Despite 
these precautions, there was some build-up of material ahead of the roller. 
The forming was stopped with a brake after two full rotations of the mandrel 
resulting in the contact patch of the roller mid-forming. 

3.1. Experimental Surface Measurement 

The workpiece and mandrel were removed from the lathe and imaged 
using a FARO Laser ScanArm0 controlled with Geomagic Studio S@ soft- 
ware. The resulting 3-D point cloud of the experimental flow formed tool- 
ing/workpiece surface was used to compare with the analytical solution. In 
an effort to assess the accuracy of the surface scan, both the mandrel and 
the contact patch were scanned together. The mandrel portion of the scan 
gave a mean of R m = 69.334 mm over 16,000 points. Defining the differ- 
ence between the scan and the contact measurement (69.331 ± 0.013 mm) of 
the mandrel radius as error, the Mean Square Error (MSE) of the scan was 
0.0366 mm. Equation |3H is the relationship used for the MSE calculation 
where n is the number of samples and 5d is the difference in distance. This 
measurement from the scan in comparison to the contact measurement made 
with a micrometer indicates that the FARO scan accuracy is on the same 
order SIS db conventional contact measurement. 



3.2. Comparison to the Analytical Solution 

The analytical surface was generated from process geometry listed above 
and was used to generate a surface of the contact interface in the form of a 



^ARO Technologies, Inc. 




(34) 



n 



2 Geomagic Inc. 
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Table 2: Mean 8d and MSE results of analytical to experimental surface with circumfer- 
ential and axial perturbation 



Condition 


Mean 8d 


MSE 


Baseline 


0.1614 


0.4017 


+0.4 mm circumferential 


0.1770 


0.4207 


-0.4 mm circumferential 


0.1648 


0.4059 


+0.4 mm axial 


0.1630 


0.4037 


-0.4 mm axial 


0.2607 


0.5106 



3-D point cloud for comparison. The two point clouds (analytical and exper- 
imental) were then overlaid. Since the origin of the experimental point cloud 
was unknown, a procedure was developed to determine the correct relative 
position and register the points of the experimental cloud. The final location 
was arrived at by axial and angular perturbations of the analytical surface 
origin. An assessment of the degree of fit of the analytical surface was found 
by calculating the nearest neighbour points on the experimental surface. For 
each point on the analytical surface, the experimental point cloud was parsed 
until a point with minimum distance was found. The distance between the 
analytical point and the nearest neighbour experimental point was defined 
as Sd- The mean of the distance between analytical and experimental points 
(Sd), and the Mean Square Error (MSE) were minimized by translating the 
analytical surface in the axial and angular directions to determine the final 
location. Table [2] shows a summary of how these parameters change by mov- 
ing the origin of the analytical surface relative to the experimental one. In 
order to compare the surfaces, an intermediate surface was linearly inter- 
polated through both the experimental and analytical point clouds and the 
distance between the surfaces (Soi) was found. 

Figure 7 shows a comparison of the experimental point cloud and the 
corresponding analytical contact patch. Due to the surface overlap, the an- 
alytical surface is depicted as being offset along the x axis by 15 mm from 
the experimental surface. At the best fit position, the analytical surface is 
qualitatively indistinguishable from the experimental surface. Figures 8 and 
9 show the 2-D plots of the analytical profile of the roller in relation to corre- 
sponding nearest neighbour experimental points on the xy plane (Fig. 8) and 
the zx plane (Fig. 9). These two plots show the proximity of the analytical 
surface to the experimental one. Outside of the area influence of the roller, 
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Figure 7: Numerical surface generated from the FARO LaserScan point data with the 
associated analytical surface offset 15 mm along the x axis. 



there is a lack of coherence with the analytical surface: this is due to the 
inherent bulging of material ahead of the roller both axially and tangentially. 
Figure 10 is a contour plot of the distance between the analytical and exper- 
imental surfaces, Soi- This plot shows that the peak distance, approximately 
0.35 mm, occurs along the axial entry profile due to the minor bulging of 
material as it encounters the roller. However, there is very little change in 
the circumferential 6 a gradient despite the overall surface variation which 
re- asserts the accuracy of the fit. 



4. Application 

The flow forming process can be thought of as a simultaneous combi- 
nation of both rolling and extrusion (or drawing) processes. An important 
application of the calculated roller/workpiece contact area is the use of the 
calculated planar projections A xy , A xz and A yz to determine the relative 
quantities of extrusion and rolling that are occurring. During flow forming, 
there is a set of running conditions that sometimes generates diametrical 
growth or defects. This is due to tangential deformation and, while very 
small, is approximated with the present analytical solution. iGur and Tirosh 
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Figure 8: Analytical profile of the roller in the best-fit position with the corresponding 
nearest neighbour experimental points on the (a) xy plane and (b) the location of these 
points on the relevant portion of the surface in 3D. 
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Figure 9: Analytical profile of the roller in the best-fit position with the corresponding 
nearest neighbour experimental points on the (a) zx plane and (b) the location of these 
points on the relevant portion of the surface in 3D. 
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Figure 10: Distance between analytical and experimental interpolated surfaces (Soi) plot- 
ted on the analytical surface 



(119821 ) proposed that the ratio of axial contact length to the circumferential 
contact length dictates the ratio of extrusion (or drawing for forward flow 
forming) to rolling that occurs during flow forming. A more accurate mea- 
sure of this extrusion/rolling ratio would be to consider the ratios of the xy 
and the yz projections of the roller /workpiece contact area. Therefore, it is 
important to see how these quantities vary with respect to the independent 
process variables. 



Using the baseline independent geometric variables presented by lXu et al 



(120011 ). a sensitivity analysis was conducted on the independent variables that 
define the three components of the surface area. The process values used by 
Xu et al.l (120011 ) were P = 0.6 mm/rev, t = 5 mm, tf = 3.5 mm, R m = 35 
mm, R r = 82 mm, R = 3 mm and an attack angle of a = 25°. Note that 
this set of variables dictates an 'AF (Section l2.ll) contact condition where the 
exit angle j3 does not participate. Through implementing One Factor At a 
Time (OFAT) analysis, whereby one variable is changed while holding others 
constant, the independent variables were changed individually between -50 
and +100% of the initial values with the exception of to and tf. The starting 
thickness was varied -25 and +100% as outside the lower range there is no 
contact. The final thickness was varied to -50%. The resultant effects on the 
contact area components were then calculated. These results are presented in 
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Figure 11: The percent change in the contact area and projections on major planes by 
varying a baseline starting thickness (t ), final thickness (tf), pitch (P), mandrel radius 
(R m ), attack angle (a), roller radius (R r ) and roller nose radius(_R). The results for R r 
and R m have A xyz , A xy and A yz overlaid. 



Figure 10. The overall effect on changing the independent process variables 
on both the area and the associated components is given in Table [3j 

The OFAT analysis technique is limited as, by definition, it does not 
allow for simultaneous changes in multiple variables. For the given geometry, 
however, this analysis does display the following important observations: 

• In terms of the largest effect on the overall contact area A X y Z ^ changing 
the material starting and final thicknesses and the pitch had the largest 
effect. This is also true for all of the area components, A xy , A xz and 
A yz . In order of precedence, the variables that had largest sensitivity 
on the overall area other than thicknesses and pitch were the radius of 
the mandrel, the attack angle, radius of the roller, with the nose radius 
having the least effect overall. 
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Table 3: Range of percentage change in contact area with percent change in independent 
process variables 



Variable 


% Range 


% Range 


in Area 


in Variable 


A 


A 


A 

^T-XZ 


A 

^-yz 


to 


125 


397.43 


389.83 


447.48 


462.14 


tf 


150 


206.30 


204.18 


214.27 


224.23 


p 


150 


73.35 


74.48 


146.31 


62.24 




150 


41.23 


41.21 


81.52 


41.34 


a 


150 


37.95 


54.65 


3.20 


50.99 


R r 


150 


20.45 


20.50 


87.07 


20.48 


R 


150 


8.60 


10.39 


3.85 


5.39 



• Varying the the roller nose radius had the least effect on the contact 
area as well as the A xy and A yz components. 

• The rolling component, A yz , followed the same trends as the overall area 
for changes in thicknesses, pitch, mandrel/roller radii and attack angle. 
This component decreased while the extrusion/drawing component and 
the overall area increased for larger roller nose radii. Furthermore, A yz , 
is more sensitive to the attack angle than the mandrel radius. 

• The largest effect on the A xy component, or the drawing/extrusion 
part of the deformation saw the same precedence of variables as for 
the overall contact area. This component showed the same response to 
variable changes as the overall area. 

• The tangential deformation component, A xz , is marginally more sen- 
sitive to the radius of the roller than the radius of the mandrel, and 
the roller nose radius has approximately the same sensitivity as the 
attack angle. This component also increased for larger values of pitch 
and mandrel radius, but decreased for larger roller radii and nose radii. 
A xz remained unaffected by changes in attack angle. 

• The overall contact area increased with increased variable values in all 
cases except for the final thickness value and the attack angle. This 
decrease was a linear for the former and non-linear for the latter. 
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• The effect of changing the starting thickness, final thickness and roller 
nose radius is a linear change for all area components while all others 
are non-linear. 

These findings are of practical importance to flow forming. If a worn roller 
is to be re-used after resurfacing, it may be necessary to modify the pitch 
in order to maintain the same forming geometry when the process was first 
commissioned. If a single set of rollers are to be used with different mandrels, 
it is also important from a process design standpoint so that the same forming 
zone geometry can be maintained. Furthermore, knowing the sensitivity of 
each of the variables on the overall contact and therefore deformation mode 
also permits easier troubleshooting of existing processes. 

5. Conclusion 

An analytical model of the roller /workpiece interface in flow forming has 
been developed such that it may predict the contact area. This model is ap- 
plicable to all tooling geometries for both forward and backward flow forming 
processes. Due to the general nature of the description of the geometry, the 
approach taken can be used for other rotary forming operations where a die 
or a roller is used to deform a cylindrical workpiece locally. This model has 
been compared to experimental data generated from physical modelling and 
shows excellent correspondence. Specifically, the analytical model was found 
to describe the experimental surface within 0.4 mm based on mean square 
error. 

An example of the application of the model has been demonstrated in 
the form of a OFAT sensitivity analysis applied to independent geometric 
variables determining tooling interaction. The independent geometric vari- 
ables examined were starting and final thicknesses (to, tf), forming pitch (P), 
mandrel radius (R m ), attack angle (a) as well as roller and roller nose radii 
(Rr, R). These variables were modified over a range of 125% for the starting 
thickness and 150% for all others. Specific findings showed that on the basis 
of a unit change in the respective variables: 

• to had f° ur times the effect on the the change in overall area and com- 
ponents 

• tf had 33% more of an effect 
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• P had 50% less of an effect, with the exception of the tangential defor- 
mation component, A xz which had 50% more of an effect 

• R m , a, and R r have less than a 27% effect 

• R caused the least change: less than 7% change in area 

The present work could be extended to study the multi- variant effects on 
the contact area to fully account for the geometric changes during compli- 
cated forming processes. However, geometric factors are not the only process 
parameters which govern the process. The main direction of future work is 
to link the geometric factors to other process factors such as workpiece ma- 
terial properties and tribological considerations to gain deeper insight into 
the overall process mechanics. 
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